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CLOSED FORM SUMMATION OF C-FINITE SEQUENCES 

CURTIS GREENE AND HERBERT S. WILF 

To David Robbins 



Abstract. We consider sums of the form 

n-l 

^ f Fi(am + bxj + c 1 )F 2 (a 2 n + b 2 j + c 2 ) . ..F k (a k n + b k j + c k ), 
3=0 

in which each {Fi(n)} is a sequence that satisfies a linear recurrence of degree D(i) < oo, with constant 
coefficients. We assume further that the a;'s and the a; + 6;'s are all nonnegative integers. We prove that 
such a sum always has a closed form, in the sense that it evaluates to a linear combination of a finite set of 
monomials in the values of the sequences {Fi(n)} with coefficients that are polynomials in n. We explicitly 
describe two different sets of monomials that will form such a linear combination, and give an algorithm for 
finding these closed forms, thereby completely automating the solution of this class of summation problems. 
We exhibit tools for determining when these explicit evaluations are unique of their type, and prove that in 
a number of interesting cases they are indeed unique. We also discuss some special features of the case of 
"indefinite summation," in which ai = a 2 = ■ ■ ■ = a k = 0. 
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1. Introduction 

In section 1.6 of the following assertion is made: 

All Fibonacci number identities such as Cassini's F n +iF n -i — F% = ( — 1)™ (and much more com- 
plicated ones), are routinely provable using Binet's formula: 

This is followed by a brief Maple program that proves Cassini's identity by substituting Binet's formula on 
the left side and showing that it then reduces to (— 1)™. Another method of proving these identities is given 
in in which it is observed that one can find the recurrence relations that are satisfied by each of the 
two sides of the identity in question, show that they are the same and that the initial values agree, and the 
identity will then be proved. 

The purpose of this note is to elaborate on these ideas by showing how to derive, instead of only to verify, 
summation identities for a certain class of sequence sums, and to show that this class of sums always has 
closed form in a certain sense, and that these closed forms can be found entirely algorithmically. Indeed, 
a Mathematica program that carries out the procedures that we develop in this paper can be downloaded 
from the web sites of the authors @|. 

We deal with the class of C-finite sequences (see QHl)- These are the sequences {F(n)} n >o that satisfy 
linear recurrences of fixed span with constant coefficients. The Fibonacci numbers, e.g., will do nicely for a 
prototype sequence of this kind. The kind of sum that we will consider first will be of the form J3J| below. 
We will say that such a sum has an F -closed form if there is a linear combination of a fixed (i.e., independent 
of n) number of monomials in values of the f's such that for all n the sum f(n) is equal to that linear 
combination. 

For example, consider the sum 

n-i 

f(n) = J2m 2 F(2n-j) 

3=0 

where the F's arc Fibonacci numbers. In Section 3.1 we will see how to use our method to show that f(n) 
can be expressed in the form 1)11(1 . which is a linear combination of five monomials in the F's. Hence the 
sum f(n) has an .F-closed form. 

More generally, we will consider functions F(n) satisfying a recurrence of minimal order D with constant 
coefficients, whose associated polynomial has roots 7*1, ri, ■ ■ ■ , rd, of multiplicities e\, e%, . . . , e<2, where J2i e i = 
D. Such a function may be expressed in the form 

d e m -l 

(1) F(n)= E X m, h n h (r m ) n , 

m=l h=Q 

where the r% are distinct and nonzero, and A m>em _i 7^ for all m. 
We will begin by evaluating sums of the form 

n-1 

(2) f(n) = F{am + hj + a) ■ ■ ■ F(a k n + b k j + c k ) 

j=o 

in which the a's, 6's, and c's are given integers. We assume further that, for all i, Oj > and a% + b, t > and 
at least one of these is positive. Later we will generalize this result to allow the F's in the summand to be 
different C-finite functions. The principal result of this paper is perhaps Theorem 1171 below, which proves 
in full generality, i.e., with arbitrary root multiplicities and with the _F"s in the summand being all different, 
the existence of closed forms, and exhibits an explicit finite basis for the solution space. 

It is elementary and well known that f(n) is C-finite, and one can readily obtain explicit expressions for 
f(n) in terms of the roots r m . Our first results show how to obtain formulae for sums f(n) of the form 
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as a polynomial in the f's, based on two different explicit sets of "target" monomials in the f's. Using the 
first target set, we obtain the following result. 

Theorem 1. The sum f(n) in has an F- closed form. It is equal to a linear combination of monomials 
in the F 's, of the form 

F((a x + 6i)n + i ± ) . . . F((a k + b k )n + i k ), < i v < d - 1 

(3) Vn,...,i fc ( n )^( a i n + h)F(a 2 n + i 2 ) ■ ■ ■ F{a k n + i k ), < i„ <d—l, 

where ifii 1 ,... t i k (n) denotes a polynomial of degree at most (3 = 1 + Afi, where A = max m {e m — 1} and 
H = \{i : a,i = 0}|. If F is rational-valued then there are solutions in which all coefficients to be determined 
are rational. 

We note that, if the roots are distinct, then (3=1 and the polynomials ipii,...,i k ( n ) are linear. The next 
theorem gives a closed form in terms of an alternate target set of monomials. 

Theorem 2. The sum f(n) in 0) can be expressed in F -closed form as a linear combination of monomials 
of the form 

F{n + h)F{n + i 2 ) . . . F(n + i P ), < i v < d- 1, 

(4) -iP ll ,... tlk (n)F(n + i 1 )F{n + i 2 )...F(n + i Q ), 0<i v <d-l, 

where P = (di + bi) + (a 2 + b 2 ) + • • • + (a k + b k ), Q = ai + a 2 + • • • + a k , and V'-ti ... i fe ip) denotes a polynomial 
of degree 7 = max{0, 1 + A(fc — J^{ai|ai > 0})}. If F is rational-valued then there are solutions in which all 
coefficients to be determined are rational. 

For example, when F(n) is the nth Fibonacci number, Theorem 2 states that any sum of the form ((21 can 
be expressed as a linear combination of monomials in F(n) and F(n + 1), with rational linear polynomial 
coefficients, where those monomials have at most two distinct degrees. Again, we note that, if the roots Tj 
are distinct, then 7 = 1 and the polynomials ipii,...,i k ( n ) are linear. 

The natural domain for these questions is the vector space V°° of complex- valued functions on {0, 1,2,...}. 
However, to obtain our expansions it is only necessary to work in the vector space V M of functions on 
{1, . . . , M}, where M is the number of unknown coefficients to be determined. More precisely, we define M 
to be equal to the number of "algebraically distinct" monomials of the form n h Y[ v F^a^, + b v )n + i„) or 
n h Y[ u F(a u n + i v ) generated by © or Here we consider two monomials to be equivalent if they differ 
by a rearrangement of factors, or by a constant factor arising from cases where a v + b v = or a v = 0. (For 
examples, see Section 3.) Then we have the following result. 

Theorem 3. Let W M C V be the vector space of complex-valued functions on {1,2,..., AT} spanned 
by the monomials in where M is the number of algebraically distinct monomials generated by 0), as 
defined in the previous paragraph. Let W°° C V°° be the vector space of functions on {0, 1, . . . } spanned by 
the same monomials. If two linear combinations of monomials of type agree in W M , then they agree in 
W°°. A similar statement holds for monomials of type ^j). 

As a consequence, we can obtain expressions of type (J3J) or Q by equating M values of f(n) to the values 
of the assumed linear combinations, and solving for the coefficients. We note that M < ((3 + 2)d k in case 
and M < d p + (7 + l)d Q in case (QJ. 

In general, F'-closcd expressions are not unique. For example, we may add terms of the form ^(F)(F(n-\- 
2) — F (n + 1) — F(n)), where ty(F ) is any polynomial in the F(an + i), to an expression involving Fibonacci 
numbers and get another valid F-closed form. However, the formats described by J3| and (0J arc highly 
restrictive, and the resulting expressions can be shown to be unique in a surprising number of cases. We 
will return to the question of uniqueness and, more generally, to the problem of computing dim(W°°), in 
Sections 4-6. 
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The structure of this paper is as follows. In Sections 2-6 we consider the case where the roots all have 
multiplicity one. In this case, both the statements and proofs of our results are considerably simpler, and 
will serve as models for the more general case to be presented later. Section 2 gives proofs of Theorems ^ 
[5] and |3 in this special case. Section 3 illustrates these results with several examples. Sections 4-6 consider 
the issue of uniqueness and dimension, again for the distinct root case. Section 7 drops the assumption 
of distinct roots, and gives a proof of Theorems Q an d in a more general form (Theorem 117(1 where the 
factors in the summand of J5J may involve different F's. Section 8 considers some issues that arise when 
ai = Oa = • • • = Ofc = 0, i.e., when the problem of computing f(n) is an indefinite summation problem. 
Section 9 contains more examples. 



2. Proofs in the case of distinct roots 

In this section, we will assume that F(n) satisfies a recurrence of minimal order d, with distinct roots, 
and hence can be expressed in the form 

d 

(5) F(n) = ]T A ro r™ 

m— 1 

with the r m distinct and the A m nonzero. Expanding the right side of above and using JSJ), we find that 

n-1 k ( d ~\ 

/w=En e«" +w+c ' • 

j=0 £=1 {m=l ) 

A typical term in the expansion of the product will look like 

(f{\ T^ r ain+b 1 j+c 1 a 2 n+b 2 j+c 2 r a k n+b k j+c k 

V u / Ji ' mi ' m 2 " " * ' m k i 

in which K is a constant, i.e., is independent of n and j, which may be different at different places in the 
exposition below. Since we are about to sum the above over j = . . . n — 1, put 



= r bl r b2 



r b k 



because this is the quantity that is raised to the jth power in the expression ©. Now there are two cases, 
namely = 1 and 0^1. 

Suppose = 1. Then the sum of our typical term (0 over j = . . . n — 1 is 

(7) Kn(rZrZ---r% k ) n - 

On the other hand, if =/= 1 then the sum of our typical term © over j = . . . n — 1 is 

(8) K { (C+ bl ■ • ■ C\ +bfc )" - (rZ • ■ ■ CJ"} • 

The next task will be to express these results in terms of various members of the sequence {F(n)} instead 
of in terms of various powers of the r^'s. To do that we write out (JSJ for d consecutive values of n, getting 

d 

F(n + i) = ]TA m C+* (i = 0,l,...,d-l) 

m— 1 
d 

- E( A ™<>™ (i = 0,l,...,d-l). 

m— 1 

We regard these as d simultaneous linear equations in the unknowns {r™, . . . ,r^}, with a coefficient matrix 
that is a nonsingular diagonal matrix times a Vandcrmondc based on distinct points, and is therefore non- 
singular. Hence for each m = 1, . . . , d, r™ a is a linear combination of F(n),F(n + 1), . . . , F(n + d — 1), with 
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coefficients that are independent of n. Thus in eqs. JJJ), ijBJ) we can replace each by a linear combina- 
tion of _F(aj7i), F(a,in + 1), . . . , F(a,in + d — 1) and we can replace each r^f*" 1- by a linear combination of 
F((a,i + bi)n), F(( ai + b t )n + 1), . . . , F(( ai + b^n + d - 1). 

After making these replacements, we see that the two possible expressions Q, (|SJ| contribute monomials 

that are all of the form ©. with the polynomials Vii i fc a U linear. This establishes the existence of 

expansions in monomials of type @, as claimed in Theorem 1. 

To prove the corresponding claim made in Theorem[21 it suffices to observe that, in the above argument, we 
could have written r^ 0i = (rj^. ) ai and replaced it by a homogeneous polynomial of degree a-i in F(n) 7 F(n + 

I) , . . . , F(n + d — 1). Similar reasoning applies to rm^' +b '\ Thus all of the resulting monomials are of type 
@. 

We continue now with the proof of Theorem |3] The arguments are identical for cases Q and (@J, so we 
will consider only case We have observed that for each i = 0, . . . , d — 1, F(n + i) is a linear combination 
of r™, . . . , rj, and conversely. Hence, in both V M and V°°, the linear span of the set 

F(ain + i 1 )F(a 2 n + i 2 ) ■ ■ ■ F(a k n + i k ), < i v < d - 1 

is equal to the linear span of the set {9™, 9%, . . . , 9^ k }, where the 9j range over all monomials of the form 

r "l r a-2 . . . „a k 
1 mi ni2 rrik ' 

Arguing similarly for the other cases, we sec that the linear span of all monomials of type J2J is equal to the 
linear span of the set of 3d k functions 

on nn rtn 

t) 1 , U 2 > • • • j o d k 

(9) 

n*l>i,mp2,...,mp% h , 

where the 9i are as defined above and the tpj range over all monomials of the form 

r ai+bl r a2+f>2 . . . „a k +b k 
mi ni2 Vtlk 

We claim that the number of distinct functions appearing in is less than or equal to M . Indeed, it is 
straightforward to check that the map F(9n + i) t-y (rf +1 )" extends to a well-defined, surjective map from 
the set of equivalence classes of monomials of type (j2J to the set of functions appearing in © • 

Now suppose that $(n) and \f(n) are linear combinations of monomials of type ©, with $(n) = \l/(n) 
for n = 1, 2, . . . , M . We know that $(n) and ^{n) can both be expressed in the form 

i j k 

for some constants c,;, dj, e k , where the sum is over distinct elements of (0 and hence there are at most M 
terms in the sum. It follows from standard results in the theory of difference equations (e.g., see [3], Chapter 

II) that <I>(n) and satisfy the same linear recurrence of order at most M with constant coefficients, 
e.g., the recurrence with characteristic polynomial Yli(t ~ 0i)Tij(t — V'j) 2 - Hence the values of $(n) and 
ty(n) are completely determined by their values for n = 1,2,..., M, and since they agree for these values, 
they must agree for all n. This completes the proof of Theorem □ 

The proof also shows that the C-finite degree of f(n) is bounded by 3d k . Sharper bounds appear in 
Corollaries El an d EOl below. 

3. Examples 

3.1. A Fibonacci sum. This work was started when a colleague asked about the sum 

n-l 

(10) /(n) = ^F(j) 2 F(2n-j), 

3=0 
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in which the F's are the Fibonacci numbers. If we refer to the general form J5J) of the question we see that 
in this case 

k = 3; d = 2; (oi, &i, Ci) = K &a, c 2 ) = (0, 1, 0); (a 3l &a, c 3 ) - (2, -1, 0). 

If we now refer to the general form yj of the answer we see that the sum f(n) is a linear combination of 
monomials 

nF(2n), F(2n), nF{2n + 1), F(2n + 1), F(nf, F(n) 2 F(n + 1), F(n)F(n + l) 2 , F(n + l) 3 . 

Hence we assume a linear combination of these monomials and equate its values to those of f(n) for n = 
0, 1, . . . , 7 to determine the constants of the linear combination. The result is that 

(11) f(n) = - (F(2n) + F(n) 2 F(n + 1) - F(n)F(n + if + F(n + if - F{2n + 1)) . 

This formula is expressed in terms of monomials of type ©. Using monomials of type Q, wc obtain the 
alternate expression 

(12) f(n) = -(2F(n)F(n + l) -2F(n f - F(n+1) 2 + F(n) 2 F(n + l) 

-F(n)F(n + lf + F{n + lf) 

In Section 5 we will show that both of these expression are unique, i.e., (|llfl is the unique -F-closed formula 
for f(n) of type © and l|12f) is the unique F-closcd formula of type (0J. 

3.2. An example involving subword avoidance. Given an alphabet of A > 2 letters, let W be some 
fixed word of three letters such that no proper suffix of W is also a proper prefix of W. For example, W = aab 
will do nicely. Let G(n) be the number of n-lcttcr words over A that do not contain W as a subword. It is 
well known, and obvious, that 

(13) G(n)=AG(n-l)-G(n-3), 

with G(0) = 1, £r(l) = A,G(2) = A 2 , so this is a C-finite sequence. It is easy to check that the roots 
of its associated polynomial equation are distinct for all A > 2. Suppose we want to evaluate the sum 
2(n) = £"=oG(i) 2 .Then 

k = 2; d = 3; (oi, b x , ci) = (a 2 ,b 2 , c 2 ) = (0, 1, 0). 

Using either Theorem 1 or Theorem 2, we see that g{n) is a linear combination of the monomials 

1, n, G(nf, G(n)G(n + 1), G(n)G(n + 2), G(n + l)G(n + 2), G(n + if, G{n + 2) 2 . 

As before, we assume a linear combination of these monomials with constants to be determined, and we 
equate the result to computed values of g(n), for n = 0, 1, . . . , 7, to solve for the constants. The end result 
is that 

(14) g{n) = ( 1 - {A - lfG(nf - 2G(n)G(n + 1) + 2G(n)G(n + 2) 

+2(A - l)G(n + l)G(n + 2) - (A - l) 2 G(n + i f - G(n + 2) 2 

if A > 2, and 

(15) g(n) = n + 2G(nf + 7G(n)G(n + 1) - 5G{n)G(n + 2) 

~5G(n + l)G(n + 2) + 3G(n + if + 2G(n + 2f 
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In the case A = 2, it is easy to show that G(n) = F(n + 3) — 1 for all n, where F(n) is the nth Fibonacci 
number. Consequently, G(n + 2) — G(n + 1) — G(n) = 1, and adding any multiple of the relation 

(G(n + 2) - G(n + 1) - G(n) - l) 2 



(16) 







to the right side of l|15|l gives another degree 2 expression of type © or Thus formula l|15fl is not unique 
within the class of formula? of type J2J) or However, in Section 5 we will show that, when A > 2, formula 
(11411 is unique within this class. When A = 2, we show that all relations are constant multiples of i|16[l . 

3.3. Fibonacci power sums. Theorems 1 and 2 imply that if the F(j)'s are the Fibonacci numbers then 
for each integer p = 1,2,... there is a formula 

n— 1 p 

f(n) = J2 F(jf - ]T A,,,, Fl /,.•'/•;/, + l)^ 1 + c p n + d p . 
j=a j=o 

Here is a brief table of values of these coefficients. 



P 







A p,2 


A P ,3 


Ap,4 


Ap,5 


Ap,6 


A p,7 




dp 


1 


























-1 





1 


-1 























i 

2 


3 
2 





3 
2 

















1 

2 





2 

25 


3 

25 


14 

25 


19 

25 











(j 

25 





7 

22 


5 

22 


15 
11 


10 
11 


15 
11 


15 
22 











7 
22 





1 

2 


5 
4 





5 
4 


1 

2 


-1 











139 
638 


763 
<,:;n 


945 
638 


350 
319 


105 

58 


357 
319 


105 
319 


777 

638 





139 
i,:;n 



The resulting expressions for f(n) turn out to be unique within the class of type @ or formulas when 
p ^ mod 4. When p is a multiple of 4 (for example, in the fourth line of the table above) the formula! 
are not unique, but are subject to a one-parameter family of relations generated by powers of the degree-4 
relation 

(F(n + l) 2 - F{n) 2 - F(n)F(n 
We will establish these facts in Section 5. 



I)) 2 



1 



3.4. Generic power sums. Consider the power sum 



f(n) = J2F(kf, 



k=0 



where F(n) solves a linear recurrence 



F(n) = AF(n - 1) + BF(n - 2) 

with initial values F(0) and F(l), where A and B are sufficiently general to insure that, if r\ and are the 
associated roots, then r\ and T2 are distinct and none of the monomials r\,r \T2, and r\ equals 1. This is 
equivalent to assuming simply that A 2 + AB 7^0, A ±(B — 1), and B^-l, 
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Using techniques introduced earlier, we can express f(n) as a linear combination of F(n) 2 , F(n)F(n + 
1), F{n + l) 2 , and 1. The solution may be computed explicitly in terms of A, B, F(Q), and F(l), and we find 
that f(n) equals 



(17) 
where 



(1 - B - A 2 {\ + B))F{nf + (2AB)F(n)F(n + 1) + (1 - B)F(n + if - K 

(A 2 - (B ~ 1) 2 )(B + 1) 



K = (1 - B - A 2 (B + 1))F(0) 2 + (2AB)F(0)F(1) + (1 - B)F(1) 2 . 

In fXZjl ; we observe a curious phenomenon: since F(n) depends on F(0) and F (1), we might expect that 
our linear equations would have led to a solution in which each of the coefficients depends on F(0) and 
F(l). However, this dependence appears only in the constant term. The next theorem demonstrates that 
such behavior is typical for power sums of C-finitc functions in which the terms in J2J) containing n are not 
present, i.e, in cases where no monomial in the roots equals 1. 

Theorem 4. Suppose that {F(n)} n >o is a C-finite sequence determined by a recurrence of order d together 
with initial conditions F(0), F(l), . . . ,F(d — 1). Suppose that the recurrence polynomial has distinct roots 
ri,...,rd, and suppose that no monomial of degree p in the Ti equals 1. Let f(n) = J2 7 j=o F(J) P > where p is 
a positive integer, and let 

(18) /(»)= E A il ,i 2 ,...,i d F(n + i 1 )F(n + i 2 )---F(n + i d ) + K 

0<H,22,---,2d<^— 1 

be the expansion of f(n) obtained according to the method given in Section 2. Then the coefficients j 3 
in JZ2J) do not depend on F(0), F(l), F(d - f ). 



Proof: Suppose that F(n) = J2 m =i ^™- r m- Define 

( Air i \ 

X(n) = 



and Y(n) = 



V ^ ) 



( F(n) \ 
F(n + 1) 



\ F(n + d- 1) / 



Y(n) = VX(n) and X(n) = V -1 Y(n) 



Then we have 
(19) 

where V is a Vandermonde matrix in the rj. It follows from (|19|l that the terms A m r^,l < m < d, can 
be expressed as linear combinations of the functions F(n + i), with coefficients that do not depend on 
F(0), F(l), . . . , F(d — 1). Using the method of Section 2, we can compute 



/(») = £ W 

3=0 



n—1 / d 

E ( E ^ 

j=0 ^ m=0 



m' m 



E( E kk ■ ■ ■ K( r ii r i2 ■ ■ ■ r i p y 

j—0 0<ii.i2;---,ip<d 

x A x \ r ii r i2 ' ' ' r i P ) n ~ 1 



(20) 



E 

0<ii,i2,...,i P <d 



" {n x r i2 ■ --r ip ) - 1 
(Krl){K 2 rl)---{K p rl) 

[ns^ ■■■r ip ) - 1 



K 
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where K is a constant. Using (|19|l . we can express all of the terms in 11201) except K as a linear combination 
of monomials in the F(n + i) with coefficients that do not depend on F(0), F(l), . . . , F(d — 1), as claimed. 
□ 

We will return to this subject in Section 8, where we prove that a more general version of Theorem 0] 
holds even when the roots rj are not distinct. 

3.5. Computational issues. We have seen in the above theorems and corollaries that we can decide the 
uniqueness of representations of certain sums in closed form if we can decide whether or not the N — 1 ) 
formally distinct monomials of degree p in the roots fi, . . . ,r<j actually are all different, when evaluated as 
complex numbers. We note here that there are various ways in which this can be done without computing 
the roots. 

For example, the elementary symmetric functions of these N monomials in the r^s are symmetric functions 
in the r,'s themselves. Since any symmetric function of the roots of a polynomial can be computed rationally 
in terms of its coefficients, the same applies to these. Once the elementary symmetric functions of the N 
monomials of degree p have been computed, the discriminant of the polynomial whose coefficients they are 
can be computed in the usual way. Thus, our condition on the roots of F can be tested without finding the 
roots, ft would be interesting to investigate in general this "hypcrdiscriminant" of degree p that is attached 
to a polynomial /, particularly with regard to how it factors when expressed in terms of the coefficients of 



4. Uniqueness and dimension: Fibonacci power sums 

In the next two sections, we investigate the uniqueness of the expansions guaranteed by Theorems ^ and 
El Motivated by Example 3.3, we first consider this question for expansions of the form 

3=0 

and, more generally, 

p 

J2 Ap,jF( n yF(n + ly-i + c p n + d p , 

3=0 

where F(n) denotes the nth Fibonacci number. In Section 5 we develop tools to help answer these questions 
for more general linear recurrences, and for other summations such as those arising in Examples 3.1 and 3.2. 
The techniques in these two sections can be viewed as refinements and extensions of the ideas introduced in 
Section 2 to prove Theorems and |3J 

Theorem 5. Let V = V°° denote the vector space of complex-valued functions on {0, 1,2,...}, and let W p 
denote the subspace o/V spanned by functions of the form F(n) l F(n + l) p ~ 4 for i = 0, . . . ,p, and let W+ + 
denote the subspace spanned by the same monomial expressions together with with the functions g{n) = n 
and h{n) = 1. Then 

(1) dim(Wp) =p+l, and 

(2) dim(W++) = ( P + 2 ' is dimsible b y 4 >' 

p [_ p + 3, otherwise. 

Corollary 6. The functions F(n) l F(n + l) p_4 , 1 < i < p are linearly independent, and the set 

{F(nYF(n + l) p -*}i<i< p U {n, 1} 



is linearly independent unless p is divisible by 4 ; in which case there is a single relation among its elements. 
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Proof: Let r\ = (1 + Vo)/2 and r 2 = (1 — \/5)/2 denote the roots of the Fibonacci recurrence polynomial. 
As noted earlier in the proof of Theorem 1, r™ and may be expressed as linear combinations of F(n) 



and F(n + 1) and vice versa. Consequently W p is the linear span of : 



'2 ' ' 



0, .. . ,p, and to prove 



statement (a) it suffices to show that these functions are linearly independent. But this follows immediately 
from the fact that the numbers r\r^ 



CP- 



are distinct, for i = 0, . . . , p. 
To prove part (b), consider the (p + 3) x (p + 3) matrix M p whose ith column is equal to the vector 

= 0, . . . ,p, and whose last two columns are the vectors (1, 1, . . . , 1) 
,p + 2). For example, when p = 2 we have 



(l,^,0f,...,0f +2 ), where 9 
and (0, 1 



/ 


1 


1 


1 


1 


\ 




r 2 
' 2 


nr 2 


r? 


1 


1 




r 4 


r\r\ 


rf 


1 


2 




„6 
' 2 


r\r\ 


4 


1 


3 


V 


' 2 


4 4 


r? 


1 


4 V 1 



Mo 



Note that detM p is the derivative at t = 1 of the (p + 3) x (p + 3) Vandermonde determinant det M p {t) 
where M p (t) is the matrix whose first p 
(l,£,i 2 ,...,£P+ 2 ). We have 



2 columns are the same as those of M p , and whose last column is 



detM 



i> 



— dot M p (t) 



d 



t=i 



n fa- *) n n c*-^) (*- r 

0<i<j<P 0<i<p 0<i<p 

n (^-^) n (!-^) 2 



0<i<j<P 



0<i<P 



It follows that det M p = only when t = 1 is a multiple root of det M p (t), i.e. 



1'2 



1 for some i. Using 



the fact that rir^ = — 1, it is easy to show that this property holds if and only if p is a multiple of 4. Thus, 
when p is not a multiple of 4, the columns of M p are linearly independent and we have dim(W+ + ) = p + 3. 

If p is a multiple of 4, then M p contains exactly two columns of Is. If one of these columns is suppressed, 
the argument just given shows that the remaining columns arc linearly independent. Hence rank(M p ) = p+2 
and dim(W+ + ) > p + 2. Since the dimension is clearly at most p + 2 in this case, the theorem is proved. □ 



5. Uniqueness and dimension: other recurrences with distinct roots 

Analogs of Theorem [S] hold for more general recurrences with distinct roots, but the exact statements 
depend on properties of the associated roots. The following theorem concerns relations among monomials 
of type and allows precise dimension computations in many cases. 

Theorem 7. Let F(n) be a solution to a linear recurrence of order d whose associated roots r±, ■ ■ ■ , rd 
are distinct, and let p and q be distinct positive integers. Let W p denote the subspace of V = V°° spanned 
by all degree p monomials of the form 

FinY 1 F(n + If 2 ■■■F{n + d-l) ld 

where i\ + i 2 + ■ ■ ■ id = P and ij > for all j . Let denote the subspace spanned by all degree q monomials 

FinY'Fin+lY 2 ■■■F(n + d-lY d and 
nFinf 1 F{n + l) 12 • • • F(n + d - lf d 
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where i\ + 12 + • • • id = q and ij > for all j. And, finally, let W++ = W p + W+ denote the subspace 
spanned by all of the above monomials. Then 

dim(W p ) = \S P \ , dim(W+) = 2\S q \ , and dim(W++) = \S P \ + 2\S q \ - \S P n S q \ 

where S p = {r^r 1 ^ ■ ■ ■ r l f \ i\ + ii + ■ • • id = p} and S q = {r^r 1 ^ ■ ■ ■ r l f \ i\ + ii + ■ • • id = q} are the sets of 
monomials in the Ti of degrees p and q, respectively, both viewed as subsets of the complex numbers. 

Corollary 8. The sets of monomials generating W p , W+ , and respectively, are linearly independent 

if and only if evaluations of formally distinct monomials in the sets S p , S q and S P U S q yield distinct complex 
numbers. 

The proof of Theorem is analogous to that given for Theorem but more careful analysis is required. First 
consider the case of W p . As noted in Section 2, each of the functions F(n), F(n + 1), . . . , F{n + d — 1) lies in 
the linear span of r", , . . . , r^, and conversely. Hence W p is spanned by the set {£>", 6% , • • • , where 

m(p, d) = and the 6j range over the m(p, d) formally distinct monomials of degree p in n, T2, . . . , r<j. 

Similar reasoning shows that W+ is spanned by the 2m(q, d) functions 

(21) ^ 2 V--,^ M) 

where the range over all formally distinct monomials of degree q in fx, . . . , and finally, is spanned 

by the m{p 1 d) + 2m(q, d) functions 

an nn rtn 
1 5 2 5***5 a m{p,d) 

(22) 

^1,^2.* ** »»^m(«,d)> 

where 6i and ipj are defined as above. TheoremHis now an immediate consequence of the following Lemma. 

Lemma 9. Let Wi,u>2, ■ ■ ■ , w m be complex numbers, and let ai, 0,2, • • • , o m be positive integers. Then the 
functions 

v?bf> l<i<m, < j < ai - 1 
are linearly independent if and only if the u)i are distinct. 

Lemma is a standard component of the classical theory of finite difference equations (e.g., JH], Chapter 
11), indeed it is the justification for the usual method of solution of such equations. It is easy to give a 
direct proof via generating functions, or, alternatively, one can give a Vandermonde-type proof based on the 
following elegant determinant formula ([2], but also see [5] for an extensive history of this formula). 

Theorem 10. Let x\, X2, ■ ■ ■ , x„ be indeterminates, and let 0,1,02, ... , a n be positive integers with = N. 

For all t, and for any integer k > 1, let 

PN (t,k) = -^(l,t,t 2 ,...,t N - 1 ) 

Let M{fli, 02, . . . , a n ) be the N x N matrix whose first ai rows are p^(xi, 0), . . . , p^(xi, a\ — 1), and whose 
next a 2 rows are Pn{%2, 0), . . . , Pn(,%2, 02 — 1), and so forth. Then 

n 

detM(ai,...,<) = n(oi-l)!!! (xj - x q ) a ^ 

where k\\\ denotes l!2! • • • k\ and OH! = 1. 
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For example, 



M(l,2,3) = 



/ 1 


•' : 1 


1 


x\ 


4 


1 ^ 


1 


■r-2 




x 2 


. 1 ■-) 







1 


2x 2 


3x 2 


Axl 




1 






x 3 


x 3 







1 


2x 3 


3x1 


4^3 




I o 





2 


6x 3 


12x1 


2Qx\ ) 



and 

detM(l,2,3) = 2{x2-x 1 ) 2 (x3-x 1 f{x3-x 2 f. 
Theorem describes relations among closed form expressions of type but the proof also yields similar 
results for expressions of type © ■ 

Corollary 11. Let F(n) be a solution to a linear recurrence of order d whose associated roots are dis- 
tinct. Let denote the space spanned by monomial functions of type |3J). Then dimWJj = 15*1 + 2\T\, 
where S is the set of all monomials of the form t^t^ 2 ■ • -t^ k and T is the set of monomials of the form 
^ai+bi_i_w2+b2 _ _ -j-Ok+bk an ^ j or eac fo ( ) i i j s one j f/jg roots ri,r2, ■ ■ ■ , J",*. The monomial functions of type 
10) are linearly independent if and only if formally distinct monomials in SUT correspond to distinct complex 
numbers. 

We omit the proof, which is analogous to that of the proof of Theorem We note that the set S U T in 
Corollarvllllis a subset of the set S p U S q appearing in Corollary |H1 and thus we obtain the following result. 

Corollary 12. Under the assumptions of Corollary if the monomial functions of type O) are linearly 
independent, then so are the monomial functions of type 0). 

We will now apply these results to some of the formulae in Sections 3.1 and 3.2. 

Corollary 13. For the Fibonacci sum f(n) appearing in Q10\). equation gives the unique F-closed 

formula of type ^) and gives the unique F-closed formula of type Q). 

Proof: By Theorem and Corollary El we need only check that, if v\ and r 2 denote the roots of the 
Fibonacci recurrence, then 

r\,r\ri,i\,r\,r\ri,rir\, and r\ 
are distinct real numbers. This is an elementary calculation. 

Corollary 14. For the sum gin) = 2J=o G{j) 2 arising in the subword avoidance problem with A = 2, 
solutions g(n) of type ^ are all given by pi us constant multiples of relation \16\) . 

Proof: The roots of the recurrence equation t 3 — 2t 2 + 1 = are ri,r2,T3, where r% and r 2 are roots of 
the Fibonacci recurrence and f3 = 1. By Theorem the dimension of the space W^q spanned by the six 
degree- 2 monomials in G(n), G(n+ 1) and G{n + 2) together with 1 and n is equal to IS2I + 2|So| — \S 2 n Sq\, 
where S2 = {rf, rf, 1, T\r<i, r\, r2}, and So = {1}- Elementary calculation shows that this dimension is equal 
to 7, hence the monomials generating W^~q are linearly independent apart from a one-parameter family of 
relations. 

Next we consider the case A = 3, as a warmup for the general case A > 2. 

Corollary 15. For the sum g(n) = 52j=o G(j) 2 arising in the subword avoidance problem with A = 3, 
formula \1$ gives the unique G-closed formula of type |3J). 

Proof: Here the recurrence equation is t 3 — 3£ 2 + l = 0, which has roots n = 1+n+r] 17 , r 2 = 1+rf , r3 = 
1 + ?/ 5 + ?] 13 , where 77 = e 27 ™/ 18 is an 18th root of unity. Again, by Thcorcm0 the dimension of the space of 
monomials W^.q is equal to |5 2 | + 2|5 | - \S 2 n Sol, where S 2 = {r\, r 2 , r§, r\r 2 , rir 3 , r 2 r 3 }, and S = {1}. 
A slightly less elementary calculation shows that the formal monomials in S2 U So are distinct, so that 
dim(W"2 ) = 8 and the monomial functions generating are linearly independent. 
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Corollary 16. For the more general power sum g(n) = X)j=o G(j) p arising in the subword avoidance 
problem, with p > and any A > 2, solutions of type are unique if and only if p ^ mod 6. 

Proof: An argument analogous to the calculation in Section 3.2 shows that formulae of type exist 
expressing g(n) as linear combinations of monomials in G(n), G(n+1), and G(n+2) of degree p, together with 
1 and n. We need to compute the dimension of W^J , which by Theorem|7| is equal to \S P \ + 215*01 ~ l-SpHSol, 
where Sq = {1} and S p is the set of all degree-p monomials in ri,r2, and r^, where ri,r 2 , and r 3 are roots 
of the recurrence equation t 3 — At 2 + 1 = 0. 

If p is not divisible by 6, the proof will be complete if we can show that formally distinct monomials 
in S p evaluate to distinct complex (actually real) numbers, and none of them equals 1. Suppose that 

r i lr 2 2r 3 3 = r i lr 2 2r 3 > : where J2 e i = Y^fi = P an( l e « fi Ior some i. Then by cancellation we obtain the 
relation r"* = rj 3 r^ k for some rearrangement of the indices, with Ui,Uj,u k > and at least one of these 
exponents positive. Using the relation r\r-2.rj, = —1, if necessary, to eliminate one of the roots, we obtain 
(after possibly reindexing), rj" 4 = ±r^ 3 with Vi, Vj > and at least one of these exponents positive. 

It is a straightforward exercise to show that the roots ri,r 2 and r% are all real, and that, if they are 
arranged in decreasing order, then r\ > 1, < r% < 1, and — 1 < 7-3 < 0. From elementary Galois theory we 
know that there exists an automorphism $ of the field K = Q(ri, r 2 , r 3 ) such that $ : n 1— > r2 ► rz 1— > n, 
i.e., it permutes the roots cyclically. Hence the equation r"* = ±.r^ holds for all three cyclic permutations 
of the roots. At least one of these equations leads to a contradiction, since |ri| > 1 and |r 2 |, |r 3 | < 1. This 
proves that formally distinct monomials are distinct, and it remains to show that none can equal 1. 

Suppose that r^r^r^ = 1, and the exponents e 4 are not all equal. Applying the identity rir 2 r 3 = — 1 
we obtain a relation of the form r^r^ 1 = ±1 for some pair of distinct i, j, with Ui, uj > and at least one 
positive. Again, this relation holds for all cyclic permutations of the indices, and consideration of absolute 
values leads to a contradiction in at least one case. Consequently, we must have e\ — e 2 = e% = e for some 
e. From the relations r\r?r?, = —1 and (rir 2 r3) e = 1 we conclude that e is even, which implies that p is a 
multiple of 6. This completes the proof that monomials in the G are linearly independent when p ^ mod 6. 
When p = 6m the relation (rir 2 r3) 2m = 1 gives relations in the G of degree 6, and so the proof of Corollary 
Hoi is complete. □ 

6. The general case of multiple roots 

In this section we show the result of dropping the assumption of distinct roots. We also consider a 
somewhat more general summation problem, viz. 

n-l 

(23) f(n) = ^ Fi(ain + b ± j + ci)F 2 (a 2 n + b 2 j + c 2 ) • • • F k (a k n + b k j + c fe ), 

3=0 

in which the factors of the summand may be different C-finitc functions. The analysis in this general case 
is similar to that in the case of distinct roots, but some additional machinery is required. The main result 
is the following, which is a generalization and also a strengthening of Theorem ^ 

Theorem 17. Let F±, F2, . . . , F k be given C-finite sequences. Suppose that, for each i = 1, ...,k, Fi(n) 
satisfies a recurrence of minimal degree D(i) whose polynomial equation has d(i) distinct roots. Denote these 
roots by r± , . . . , rjj2j , and let , . . . , be their respective multiplicities, so that D{i) = . Finally, 

let A(i) = maxi<j<<j(i)(e^ — 1). Then the sum f(n), of can be expressed as a linear combination of 

the monomials 

(24) F 1 ((a 1 + b 1 )n + i 1 )...F k ((a k + b k )n + i k ) (0 < «„ < D{v) - 1; 1 < v < k), and 

(25) i>h,-,i k (n)Ji(am + ii) . . . F k (a k n + i k ) (0 < i„ < D(v) -l\l<v<k), 
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in which ipi lt ...,i k (n) is a polynomial whose degree is bounded above by 

• zero, i.e., the factor ip can be omitted, if for all sequences (m 1 ,m 2 , . . . ,m k ) with 1 < m, < d(i), we 
have Yli =1 {rm i ) bi ^ 1, and 

• 1 + : a i = 0} otherwise. 

Proof: We have, for the sum f{n) of (|23|) . 

n-i fc d(^) e»-l 
(26) f(n) = J2IiH E O a ' n + fe « + 

j=0 £=1 m=l h=0 



where the coefficients X^h arG defined by the form of the i^'s, namely 

m=l h=0 

If we expand the product and the two inner sums in i|2tj|) , we find that a typical term is of the form 

K(a in + b,j + Cl ) hl (rW)° in+ ^ +Cl . . . (a k n + b k j + c k ) h * ( r (J>)«*«+6«+<* , 
where < hi < em, — 1 for 1 < £ < fc. If we write 

(a,n + &jj + a) hi = (a,i(n - j) + (a l + h t )j + c l ) hi 
and further expand each of these factors, we find that our typical term now can be expressed as 

(27) Kf(n - jY ((r«) 61 . . . (r«) h *) J ((r«r ■ ■ ■ (^) a& )" , 
in which 

(28) q<J2{ h i :a i + b i^0}, r < ^{hi : a* £ 0}, and q + r < ^h t < ]^(e$. - 1). 

i i 

At this point we need the following result. 
Lemma 18. We have 

V 7 -f„ - rtj = \ Pa{n)xn + l/ - T ^ 1: 

^ ' \^+6 + i(n), z/x=l, 

where P s (n) denotes a generic polynomial of degree s, whose coefficients may depend on x. 
Proof. Suppose that x ^ 1. If b = 0, we have 



= I a: 



dx / ' V <ix / V x — 1 

i=o x 7 j=o 



which is of the form stated when 6 = 0. For b > we have 

n— 1 y , \ a n—1 /W\ a n /^/\ a n 

gj-O-flV = (4) gc-^-^s) g*- J -(*|) 



"i ) x " ( (Pb(n) + nb)x ~ n + K ^ = ( x i) ( Pfe(n) +nb+ Kx ^ 



which is evidently of the desired form. The case x = 1 is elementary. □ 
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If we sum the typical term (|27JI over j = 0, .., n — 1 and use the Lemma, we find that the overall sum f(n) 
is a sum of expressions of the form 



n-l 
3=0 



f{n - JY (V«) hl . . . (r^lpy ((r«) ai . . . (r«)°f 



(29) 



K ((r$r ■ ■ ■ (r«)-) n x {™ " " ' + if 9 * 1; 

V " 7 tP 9+r+ i(n), if 9 = 1. 

•g(n) ((r, ( r !!) Ql+bl . . . (r£j>)°* +l *) n + #P r (n) ((rl^ . . . (r^) a *)" , if # 1; 
(n)((r«r ...(r^)^)", if = 1, 





where = (rm\) bl ■ ■ ■ (rm k ) bk , q and r satisfy the bounds given in l)28|l. and P s (n) denotes a generic 
polynomial of degree s. 

Considering each of the three terms appearing in the last member of l|29|) , we first have 

^W((^) ai+6l -(r(lr +tl )" = KP q (n) J] (r«)^+ b ')" 

{i: ai +6i^0} 

n (^) (ai+bi)n , 

say. Since q < Y^ii^i '■ a i + 7^ 0} — Si{ e mi — 1 : o» + 6 2 ; ^ 0}, each exponent j in the range < j < q can 
be written (in many ways) as j = j\ + ja + • • • + ife, where < j'j < — 1 and j'j = if a, + hi = 0. Hence 
the last member above may be expressed as 

(30) "*(»$) n(o,+6i) (0<ji<eW). 

3=0 {i:ai+bijtO} 

Now we observe that the solution space of the recurrence satisfied by Fi(n) has dimension D(i), and that 
the D(i) shifted sequences {Fi(n)}, {Fi{n + 1)}, . . . , {Fi(n + D(i) — 1)} are linearly independent, since in the 
contrary case the function F{ would satisfy a recurrence of degree < D(i). Consequently these D{i) sequences 
are a basis for the solution space, and therefore each of the functions r\? (rm) n ', j = 0, 1, ... , e™ — 1, m = 
1, . . . , d(i) can be written as a linear combination of -Fj(n), Fi(n + 1), . . . , Fi(n + D(i) — 1). 

Thus we return to our general term i|3U|) and we replace each of the monomials of the form n J Vm by 
such a linear combination of functions of the form F((a + b)n), F((a + b)n + 1), . . . , F((a + b)n + D(i) — 1), 
and expand everything again. Now our general term is of the form 1)240 in the statement of the theorem, 
which concludes the treatment of the first term in the final display of equation l|29|) . The second term, 



may be handled similarly. In the third case, which occurs when = 1, powers of n in P q+r+ i{n) can be 
redistributed as in (|30f) provided that each corresponding a, is nonzero. If all Sj / 0, a term of degree 
at most 1 remains, since q + r < e™, — 1; more generally, the residual exponent is at most equal to 
1 + £JA(i) : ^ = 0}. □ 

By keeping track of the number of terms of the form n? Y[ rf™ and n> Y[ ?"m + ^™ being summed in (|29|l we 
can give a bound on the C-finite degree of /(n), generalizing Theorem |3| to the case of multiple roots and 
different F, . 
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Corollary 19. Suppose that F\, F2, ■ ■ . , Fk are as defined in Theorem Let M denote the number of 
algebraically distinct monomials generated by \24\j and l&5|) . as defined as in the paragraph preceding Theorem 
Then the sum f(n) in is C-finite, of degree at most M. The coefficients of the monomials expressing 
f(n) as a linear combination of those monomials can be found by solving equations involving at most M 
values of f(n). If a solution is valid for the first M values of n, then it is valid for all values of n. 

Proof: Assume first that no expression of the form = ULi( r ™ e ) be equals 1. Define Q = {£ \ a e + bi ^ 0} 
and R = {£ \ ae 7^ 0}. The proof of Theorem 1171 shows that f(n) is in the linear span of the set UqU Ur, 
where Uq is the set of functions of the form 

(31) n^( r ^) (a " +Mn 

£(EQ 

where 1 < me < d(£) and < ji < e$ t — 1, and Ur is the set of functions of the form 

(32) U nn ( r ^r n 

where 1 < mi < d(£) and < ji < em e — 1. We claim that \Uq UUr\ < M. The argument is similar to the 
one used in Section 2 to prove Theorem^m the distinct root case: if 6 7^ 0, we define a map Fe(6n + i) i— » 
ni(rm) 6n i where i 1— > (j,m) is some enumeration of the D(£) pairs with 1 < m < d(£),0 < j < em, — 1. 
This map extends to a well-defined surjective map from the set of equivalence classes of monomials of type 
(12411 and (|25() to the set Uq U Ur. The maximum C-fmitc degree of any function in (Uq U Ur) is equal to 
the dimension of that space, and hence the C-finite degree of f(n) is less than or equal to M. 

Next suppose that there exist expressions of the form O = rifcil 7 ™*)^ equal to 1. The proof of Theorem 
1171 shows that f(n) is in the linear span of the set Uq U Wr, where Uq is as defined above, and Wr is the 
set of functions of the form 

(33) n h Y[(r£lr n 

eeR 

where 1 < me < d(£) and < h < H , where 

(34) ^ = i + E(^-D+E( e 5-i) 

l£R l^R 

and, for a given sequence mi with I € R, the m\ are chosen to maximize H over all sequences (mi, . . . , m^) 
such that UieR(^l) bt Y[ tm {r%) be = 1- 

As in the first case, the map Fi(0n + i) ^ n^(r l m ) 9n : 9 7^ 0, extends to a well-defined surjective map from 
the set of equivalence classes of type l|24() and (12,") I) to the set Uq U Wr, proving that \Uq U Wr\ < M. Hence 
the C-finite degree of f{n) is at most M in all cases, and the remaining assertions follow immmediately. □ 

It is possible to give another bound on the C-finite degree that is sometimes sharper than the one in 
Corollary HI 

Corollary 20. Suppose that F\, F2, . . . , F^ are as in Theorem \ll\ where each Fi is C-finite of degree D(i), 
with d(i) distinct roots. Then the C-finite degree of f(n) is bounded by 

where Q = {£ \ ai + be ^= 0} and R = {£ \ ai ^ 0}. This last expression is in turn bounded by 

(36) l[D(£) + l[D(e). 

eeo teR 
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In all of these expressions, the empty product is taken to be equal to 1 . 

Proof: Again, we first consider the case where no product of the form O = IIf=i( r n2) b * equals 1. To 
establish the bound in note that 

(37) \U Q \< ]T fl + E^- 1 )) = (n^)) + (E( D W II d(j)-l[d(£) 

n*>) (-Ed- 



eeQ 



where the first summation is over all sequences {mi)g^Q, satisfying 1 < mg < d(l) for £ £ Q. The second 
term in 1(35(1 is handled similarly, and the result follows. To obtain the bound in 1)36(1 . it suffices to show that 



K eeQ I \ eeQ 



To see this, divide both sides by rLeQ d(£), and write D{t)/d{l) = 1 + tt, where ti > 0. It remains to show 
that 1 + X^£eQ tl — WeeQ^ + ^) f° r nonnegative tg, which is obvious. This completes the proof of Corollary 
1191 in the first case. 

Next suppose that there are products of the form = Il<>=i( r ™f ) &f ec l ua l to 1. Let Wr be defined as in 
the proof of Corollary 1 191 We will show that \Uq U Wr\ is bounded by 1)35(1 . Arguing as in the first case, we 
obtain that \Uq\ is less than or equal to the first summand in 1(35(1 . Before computing the Wr contribution, 
it is convenient to define 

W° R = {n h Y[(r£lr n \H <h<H} 

leR 

where H = 1 + X^ e fl(4£ ~ !)■ We claim that 

\u Q nW R \ > \W%\. 

To see this, suppose that n h Y[teR( r< me) aen S W R . Since h > Hq, there must exist indices m| with £ £ R 
such that Yl( e R(rm e ) bt ]J e ^i?( r m* ) he = 1- Assume that these have been chosen so that the right hand side 
of l(34|) is maximized. Then 

eeR eeR t&R 

is an element of Uq H Wr, and it is clear that this map is injective. It follows that 

\u q uWr\ = \u Q \ + \w R \-\u Q nW R \ 

< \U Q \ + \W R \-\W R \ 

= \Uq\ + \Wr-W° r \. 

Finally, we have 

\w r -w°r\< (i + E^- 1 ))' 

(m e )e €R eeR 

which is less than or equal to the second summand in ((35(1 . by the argument presented in the first case. The 
bound (fTTfT|> follows as before, and the proof is complete. □ 

We remark that the second case above can also be derived from the first case by a continuity argument. 
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7. A STRIKING PROPERTY OF INDEFINITE SUMMATION 

When a\ = ai — ■ ■ ■ — a k = in l|23(l wc arc doing indefinite summation, i.e., the problem is equivalent 
to finding a function S(n) (called an indefinite sum function) such that 

(38) S(n) - S(n - 1) = F l {b 1 n + cj) ■ • ■ F k (b k n + c k ) 

In this section we will show that, if no products of the form rii=i( r »™i) bi are equal to 1, there exists an 
indefinite sum function S(n) expressible (formally) as a linear combination of monomials in the Fi with 
coefficients that arc independent of the initial conditions satisfied by the various -Fj. 

Theorem 0] proved this result (which wc will call the independence property) for power sums J^ - F(j) v , 
assuming that the roots n are distinct, and it is easy to extend that proof to the case of summands of the 
form F\{b\n + c±) ■ ■ ■ Fk(b k n + c k ), with different F^, as long as the roots are distinct. Among other things, 
this section extends Thcorcm0]to the multiple root case, where the proof turns out to be considerably more 
difficult. 

Similar questions have been considered in jjjj, where the author describes a method (different from ours) 
that finds indefinite sum functions in some but not all cases 1 , but docs not delineate these cases with a 
theorem. In fact, the method developed in [3] assumes the independence property; the author asserts (in 
somewhat vague terms, and without proof) that it holds if no products rii( r mi) bi are equal to 1. 

Theorem [^] (below) establishes this result, and in full generality. It dispenses with the assumption 
of distinct roots, and makes appropriate modifications in the case where products of the form rii( r mi) bi 
are equal to 1. The first step in the proof is to identify a collection of functions analogous to the terms 
■V™ r m> 1 — m — d that appear in the proof of Theorem 0] which can be expressed as linear combinations of 
the F(n + i) with coefficients that do not depend on the initial conditions. 

Lemma 21. Suppose that 

d e m -l 

F ^= E E X m, h n h (r m ) n , 

m=l h=0 

where the r m are distinct and nonzero, and \ m ^ m -\ ^ for all m. For m = 1, . . . , d and h = 0, . . . , e m — 1, 
define 

®T(n) = E QVi^-Vm)" 
i=h ^ ' 

Let D = ^2 m e m . Then each of the functions $}"(«.) may be expressed as a linear combination of the 
functions F(n), F(n + 1), . . . , F(n + D — 1), with coefficients that do not depend on the values of X m ,h, i-e., 
those coefficients do not depend on the initial conditions of F. 

Proof: Let F be the column vector of length D whose ith element is equal to F(n + i), < i < D — 1. 
Let R and $ be column vectors of length D, each indexed by (m, h), 1 < m < d, < h < e m — 1, where the 
(m, h)th element of R is equal to n h (r m ) n and the (m, h)th element of $ is equal to < i ) e™ l _i_/ l ("-)- Then we 
have 

F = MR 

where M is a matrix with rows indexed by 0, . . . , D — 1 and columns indexed by (m, h), whose entry in row 
t and column (to, h) is equal to 

(r m y E ( l h ) Xm ^~ h 

i—h ^ ^ 



The authors thank the referee for bringing this reference to their attention, and also for several other helpful suggestions 
and remarks. 
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A 



m,e m — 1 — {h — h') 



We have the factorization 

M M A 

where Mo is a matrix indexed as in M, whose entry in row t and column (m, h) is 

(r m ) t i e '"- 1 -'\ 

and A is a matrix with both rows and columns indexed by (to, h), whose entry is row (m,h) and column 
(to', h') is equal to zero if to + 1 to' or h! > h, and is otherwise equal to 

'e m — l — (h— h 1 ) 
h' 

One can also check that 

* = AR. 

We note that Mo does not depend on the values of \ m .h- Furthermore, after factoring out various powers 
of the roots r m , it is column-equivalent to the transpose of the matrix whose determinant was computed in 
Theorem If 01 Hence, since we are assuming the r m to be distinct and nonzero, Mo is nonsingular. Thus we 
can write 

Mo^'F = AR = 

and the lemma follows. □. 

As an illustration of the previous lemma and its proof, consider the function 

F(n) = ar" + bur" + cn 2 r r { + dr% + enr% 

Then 



$ (n) 


= ar n - 


f- bn r™ + en 2 r™ 




= 6r™- 


- 2cn r™ 




= cr™ 






= dr%- 


1- en r™ 


$f(n) 


= er%. 





We have 



MR 



/ 

(a + b + c)ri 
(a + 26 + 4c)r 
(a + 36 
V (a + 46 



1 

9c)rf 
16c)rf 



(6 

(6- 

(b 

(b 



2c)n 

Ac)r\ 
6c)r? 
8c)rf 



2 

erf 

cr\ 
,,1 



d 

(d + e)r 2 
(d + 2e)rl 
(d + 3e)r| 



er 2 



cr\ (d + Ae)r\ er\ J 



n r{ 

1 a 



and 



Also 



§i(n) 
§5(n) 
$?(n) 

V J 



( c 
b 
a 


V 





2c 
6 





\ 







V 



/ F(n) \ 
F(n + 1) 
F(n + 2) 
F(n + 3) 
V F(n + 4) J 





/ 








1 





1 


\ 


/ 


c 











\ 






n 


ri 


ri 


'■'2 


'2 






6 


2c 











M = M A = 




Ar 2 


2r\ 


1 


2r 2 


2 

'2 






« 


6 


c 












9r\ 


ir\ 


'"1 


3r| 


,,3 
' 2 















e 









16rf 


Ar\ 


ri 


Art 


r 4 
' 2 


/ 


V 













e / 



AR 



Thus the relation <fr = Mo 1 F expresses each $™(n) as a linear combination of F(n), . . . , F(n + 4), with 
coefficients that do not depend on a, 6, c, d, and e. 
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To prove our main result, we also need the following Lemma. 
Lemma 22. If p > 0, define S p (n,x) = J2]=o 3 V xJ ■ V x ^ ^> then 

S p (n,x) = ^W^lEfl^f 1 )) 

k=0 ^ ' 

= B p (x)-x n ((n + Q)P ) 

Q k ->B k (x) 

where Bk(x) = Ak(x)/(1 — x) k+1 , and Ak{x) denotes the Eulerian polynomial of degree k. Here, the notation 
Q k — ► Bk(x) means "replace Q k by Bk(x) throughout". If x = 1, then S p (n,x) is a (well-known) polynomial 
of degree p + 1 . 



oo 



Proof: The identity 

E^ P ~ 

3=0 K > 

is classical (see e.g. PP), and immediately implies the relation xB' p {x) = B p +x(x), as a formal power series 
identity. The lemma now follows easily by induction. □ 

The next Lemma provides an explicit form for our indefinite sums, expressing the result in terms of the 
functions <&™{n) defined in Lemma I2T1 

Lemma 23. Suppose that for each i = 1, . . . , k, 



m=l h=0 

with A*' m 7^ for all i and m. For 1 < i < k, let Si m = ef — 1. Then 

vn,ern —1 

n-1 

(39) 52F 1 V)F 2 (j)---F h V) = *(n) + 5(0)-S(n) 

3=0 

where ty(n) is a polynomial, and 

(40) S(n) = Yl ' 'E ' * B s(r£l--r%l){ E ^ (») ' • • (»)} ■ 

r w „. r w ^ x s=0 Ui+— +**=« J 

J/ no product of the form rll] • • • rml equals 1, then ty(n) = 0; otherwise ty(n) has degree at most 1 + A(i), 
where A(i) = max 1<:) < d ( i )(ej I ' ) — 1), as defined in Theorem \Vl\ In formula \4<Jj) , B s {x) = A s (x)/(1 — x) s+1 , 
as defined in Lemma W^ and for 1= 1, . . . , k, 

as defined in Lemma Wl\ 

Note that, in expression 1401 . the coefficients B s (r£j . . . rml) depend only on the roots rf, and hence by 
Lemma 12 II S(n) is independent of the initial conditions of the F^. 

Proof: To simplify notation, we will first consider the case where the Fi each have a single root, i.e. 

F(n) = (A^ ) +A»n + ... + Ag ) n 4 «)(rW)" 
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where rW is a root of multiplicity e^ 1 ' = 8i + 1. We have 

n — 1 n— 1 <5i 

E*W)-*W) = ECE^W^V-CE^W*')') 

.7=0 j=0 i=0 i=0 
n-1 <5i H h«5 fc 

= E E E ^■■■t ) f(r {1) --r^y 

j=0 s=0 iiH Mk=s 

<5x H h<5 fc n— 1 

-EE ^..-^SArW...^. 

s=0 iiH Mfe=s 3=0 

Now for any integer p > 0, define 



(41) 



H p (n,x) = x n J2[ P k \n k B p _ k (x), 



fe=0 



where Bk{x) is as defined in Lemma l2"2l Thus, by that Lemma, we have 

n-l 

J2fx j = H p (0,x)~H p (n,x)=x n (n + Q) p 
j=o 

for any x ^ 1. Continuing with the above calculation, we obtain 

n-l S!+-+S k 

E*tfj)-*«j) = E E Ag'.-Ai? (jy.(0 1 rW-rW)-H.(n J rW...r< 

j=0 s=0 iiH hifc=s 

provided that r^ 1 ' • • • A kS> ^= 1 ; otherwise the sum is equal to ^(n), a polynomial of degree at most equal to 
f + J2i Si = 1 + J2i( el — !)■ ^ n tne former case, write 

S(n) = ^E^ E f...A!fff ! (n ) r«...rW) 



S — 2lH hi'A;— S 



so that 
(42) 



js(0) - S(n) if • • -r( fc ) ^ 1, and 
l*(n) otherwise. 



n-l 
3=0 

Using the alternate form in (|41|l . we can write 

S(n) = 5l J2 5k E ^•■•AS ) (r-«---rW)> + Q) 

s=0 ii H Kfc=s 

E A?> (n + QV) (rW )" • • • ( E Af 5 (n + Qf) (r^ ) 



»B<,(rW— r<*0) 



i=0 
Si 



s _»B s ( r (i)... r (fe)) 



u— z— u 
<5i 



i=0 «=0 



E^E 

i=0 ti=0 

EE ( l )^ ] n- u (r^rQ u ) ■■■ ( EE ( ; )A| fc) n^(r«)"g" 



u=0 i=?z 



3 -»S s (r< 1 )---r( fc )) 



'->B,(rW.., r (*)) 



u=0 



u=0 
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Combining this last expression with l)42|l yields the statement of Lemma l23l in the case where each recurrence 
has a single root. For the general case (i.e., when there are multiple roots for each Fi), one can collect terms 
in the expansion of F\ (j)F% (j) ■ ■ ■ F k (J) corresponding to each choice r mi , r m2 , . . . , r mk of a sequence of roots 
from each Fi, and an expression for K(n) of the form l|39|l results, with S(n) as in (|40|l . This completes the 
proof of Lemma [33] □. 

Combining Lemmas 1211 and 1231 we obtain the following Theorem, which is the main result of this section. 
Theorem 24. If F\ (n) , F 2 (n) , . . . ,F k (n) satisfy the hypotheses of Lemma \2!A then the sum 

n-1 

Fi(hj + c 1 )F 2 {b 2 j + c 2 ) ■ • • F k (b k j + c k ) 

j=o 

may be expressed as ty(n)~S(n), where ^(n) is a polynomial and S(n) is a linear combination of monomials 
of the form 

Fi(6m + h)F 2 (b 2 n + i 2 ) ■ ■ ■ F(b k n + i k ) (0<i v < D{v) -l;l<v<k) 
with D(v) — ^2 m em\ such that the coefficients of that linear combination are independent of the initial 

conditions of the Fi. If no product (rm\) bl ■ ■ ■ (fml) bk of the associated roots equals 1, then ty(n) is a 
constant; otherwise it has degree at most 1 + £^ In general, ^(n) will depend on the initial conditions 

of the Fi . 

8. More examples 
We will give some examples illustrating the results in the previous two sections. 

8.1. A mixed convolution. Let F(n) denote the nth Fibonacci number, and let G(n) be defined by 
the subword-avoiding recurrence l|13|) with A — 3, in other words G(0) = 1,G(1) = 3, G(2) = 9, and 
G(n) = 3G(n — 1) — G(n — 3) for n > 2. Then we have the following identity: 

$3 J F(j)G{n - j) = 18G(n + 1) - (9G(n) + 5G(n + 2) + 3F(n) + nF(n) + nF(n + 1)) 

The target monomials in this case are 

F{n), nF{n), F{n + 1), nF(n + 1), G(n), G(n + 1), G(n + 2) 

and the (unique) solution is obtained by solving a system of 7 equations in 7 unknowns. Here we are applying 
Theorem 1171 in the case where no product of the form J^r^ equals 1, with Fi(n) — n,F 2 {n) = F(n), and 
F 3 (n) = G(n). The sum is C-finite of degree 7, and this degree achieves the bounds given in Corollaries IT§1 
and EH 

8.2. The independence property. We will give two examples of indefinite summations illustrating the 
phenomena described in Theorems HI and 1241 Consider the sum Y^j=o F(j) 3 where F(n) satisfies the Fi- 
bonacci recurrence with initial conditions F(0) = p, F(l) — 1. The target monomials are 

1, F(n) 3 , F{n) 2 F(n + 1), F(n)F(n + l) 2 , F(n + l) 3 

and we obtain the identity 

X>0f = i(l-3 P + 3 P 3 )-^F(n) 3 + ^(n)F(n + l) 2 -^(n + l) 3 , 

which is of the form 5(0) — S(n) in the notation of Theorem 1241 Since the associated roots are distinct and 
no product equals 1, this case is covered by Theorem 0] 
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Next consider the sum Xw=o -^(i) 4 s f° r which we now have products of roots equal to 1. The target 
monomials are 

1, nF(n) 4 , F(n) 3 F(n + l), F(n) 2 F(n + l) 2 , F(n)F(n + l) 3 , F(n + l) 4 
and we obtain the general solution 

i(A + A x n + 52F(n) 3 F(n + 1) - 22F(n) 2 F(n + if - 3QF(n)F(n + l) 3 + 19F(n + l) 4 ) 
+K{(-1 + p + p 2 ) 2 - {F(n + l) 2 - F(n) 2 - F{n)F(n + l)) 2 

where 

A = -19 + 36p+22p 2 - 52p 3 , A 1 = 6(-l + p + p 2 ) 2 , 

and X is an arbitrary constant. If p 2 +p— 1^0, the polynomial term \l/(n) in Theorem 1241 will always have 
degree 1. 

8.3. Partial summation of series. Consider the sum 

n-1 

where F(n) is the nth Fibonacci number and x is an indeterminate. The summand is a product of two 
C-fmite sequences, one of degree two and the other of degree one. Following TheoremEl we construct a list 
of target of target monomials 1, F(n)x n , and F{n + l)x n , and from these we obtain the identity 

£ FfcV = J a - x n X F(n) + X 2 F(n + 1) 

z — ' 1 — x — x z V 1 — x ~ x A \ — x — X A 

This identity quantifies the remainder term in the Fibonacci generating function (an equivalent result appears 
as problem 1.2.8.21 in [S]). Our approach can be easily extended; for example, using 1, F(n) 2 x n , F(n)F(n + 
l)x n , and F(n + l) 2 x n as target monomials and solving four equations in four unknowns, we obtain the 
partial summation formula 

< 43 » E^) , -'- 1 -J-'J + «. -^w 

3=0 

where 

(44) i?„(x) 



(1 - 2x - x 2 )F(n) 2 + 2x 2 F(n)F{n + 1) + x(l - x)F(n + l) 2 



The first term in (j43(l is the full generating function for squares of Fibonacci numbers. A formula for the 
full generating function for all powers p appears in [Hj (see also 0, problem 1.2.8.30). 

We note that, to obtain (|43(1 and 1)44(1 by this method, it was only necessary to know the first four values 
of the sum, and also that F satisfies some 2-term recurrence with constant coefficients. 

8.4. The degree of the polynomial multiplier. In Theorem 1 171 we gave a set of monomials in terms of 
which the sum can be expressed, in the general case of repeated roots. In those monomials a polynomial 
factor Vii.-.ifc appears, and the degree of that polynomial was found to be at most 1 + : a * = 0}- We 

remark here that the C-fmite function F(n) = n p , for positive integer p, shows that this upper bound can 
be achieved. For here we have 

k = 1, d(l) = 1, r{ 1] = 1, e[ 1} =p+l, A(l) = p, D(l) =p+l, (oi, h, a) = (0, 1, 0). 

The monomials in the list (|24|l are all of degree p, and those in the list 1(25(1 are of degree equal to the degree 
of ipi(n). The maximum allowable degree of the latter is 1 + A(l) = p + 1, and in this case ip is of degree 
p + 1 since the sum obviously is. 
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8.5. An example from the theory of partitions. Our algorithm for summation of C-finite sequences 
can sometimes have by-products that are more interesting than the particular problem being solved. A 
small example of this is given here. Suppose p^{n) is the number of partitions of n into < 5 parts. Then 

^ )1>0 p 5 (n)z n = ((1 - z)(l — z 2 ){\ — z 3 )(l — z 4 )(l — z 5 )) , so p 5 is C-finite of degree < 15. We asked our 
Mathematica program 0] to find f(n) = Y^o<j< n -i Psti)- I n addition to giving the answer, a number of the 
arbitrary constants that are used to form linear combinations with target monomials were left unassigned, 
and since the coefficient of every such unassigned constant must of course vanish, one has found an identity 
On this occasion we chose one symmetrical looking such identity from the output, namely 

p 5 (n - 10) + 4p 5 (n - 9) + 9p 5 (n - 8) + 15p 5 (n - 7) + 20p 5 (n - 6) + 22p 5 {n - 5) + 20p 5 (n - 4) 

+ 15p 5 (n - 3) + 9p 5 (n - 2) + 4j> 5 (n - 1) +p 5 {n) = 

The coefficients on the left side are recognized as the numbers of permutations of 5 letters that have k 
inversions. From this one might suspect that we have generally 

(45) J2b(k,j)p k (n-j) = f 

j ^ 

where b(k,j) is the number of fc-pcrmutations that have exactly j inversions and pfe(m) is the number of 
partitions of m into parts < k. A proof of this identity by generating functions is quite trivial. Here is a 
bijective proof, that is, a bijection between pairs consisting of a permutation of k letters with j inversions 
and a partition of n — j into < k parts, on the one hand, and one of the compositions of n into k 

nonncgative parts, on the other. Take such a composition X of n into k parts. Perform a "modified bubble 
sort," whereby whenever one sees an adjacent pair x y with x < y, it is replaced by (y— l)x. Keep doing this 
until there are no adjacent pairs x < y, i.e., until a partition A (perhaps followed by 0's) is obtained. Call the 
resulting permutation of positions a. Then A and a are uniquely determined by X, and the correspondence 
is bijective. □ 

We do not claim novelty for this result or its proof, but offer it only as an example of the usefulness that 
our algorithms can have in the discovery process. 
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